Self-consistent local-equilibrium model for density profile and distribution of 
dissipative currents in a Hall bar under strong magnetic fields 
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Recent spatially resolved measurements of the electrostatic-potential variation across a Hall bar 
in strong magnetic fields, which revealed a clear correlation between current-carrying strips and 
incompressible strips expected near the edges of the Hall bar, cannot be understood on the basis 
of existing equilibrium theories. To explain these experiments, we generalize the Thomas-Fermi- 
Poisson approach for the self-consistent calculation of electrostatic potential and electron density 
in total thermal equilibrium to a local equilibrium theory that allows to treat finite gradients of the 
electrochemical potential as driving forces of currents in the presence of dissipation. A conventional 
conductivity model with small values of the longitudinal conductivity for integer values of the (local) 
Landau-level filling factor shows that, in apparent agreement with experiment, the current density 
is localized near incompressible strips, whose location and width in turn depend on the applied 
current . 



I. INTRODUCTION 

The question, where the current flows in a Hall bar 
under quantum-Hall-effect (QHE) conditions, has been 
investigated by many authors, and many controversial 
answers have been given. Part of this controversy arises 
from the fact that, apart from the total electric cur- 
rent density, different partial current densities can be de- 
fined which integrate to the same total current. 1 Even 
in the thermodynamic equilibrium state with vanish- 
ing total current, quantum calculations yield finite cur- 
rent densities which are related to the density variation 
and the energy dispersion of the two-dimensional elec- 
tron gas (2DEG) near the sample edges. Such current 
distributions have been calculated using different ap- 
proaches, from simple Hartree-type approximations 2 to 
sophisticated treatments of exchange and correlation ef- 
fects within a current-density functional theory 3,4 . Un- 
fortunately, there is little experimental information about 
the current-density distribution in samples with zero or 
small total current. Experiments with high currents, of 
the order of the critical current for the breakdown of 
the QHE, have shown that the current distribution de- 
pends strongly on the mobility of the sample 5 , but also 
on screening effects caused by nearby metallic gates 6 or 
mesoscopic inhomogeneities like arrays of antidote ' . 

Our present work is motivated by a recent series of 
experiments 9 " 11 , in which a scanning force microscope 12 
is used to measure the Hall-potential distribution across 
a Hall bar under QHE conditions. The characteristics 
of the potential distribution are found to change drasti- 
cally with the magnetic field applied perpendicularly to 
the sample, i.e., with the Landau level filling factor v. 
While for v values far away from integers the potential 
varies linearly across the sample ( "type I" behavior) and 
for integer and slightly lower values a non-linear poten- 
tial drop in a broad region in the middle of the sample 
is observed ("type II"), for slightly larger than integer 



v values the potential is flat in the center region and 
drops across two strips that move with increasing v to- 
wards the sample edges ( "type III" ) . Position and width 
of these strips coincide with those of the incompressible 
strips 10 that are expected to form in the sample due to 
the non- linear screening properties of the 2DEG in strong 
magnetic fields 13 " 15 . 

The fact, that the Hall voltage drops across the incom- 
pressible strips, indicates that the current flows prefer- 
ably along these strips. 2 Such a conjecture can be found 
in an early paper by A. M. Chang 16 . Self-consistent equi- 
librium calculations, which imposed a dissipationless Hall 
current as a thermodynamic boundary condition on the 
2DEG in a Hall bar, 17 ' 18 could however not confirm this 
conjecture. On the contrary, the current-density profile 
was found to extend over the whole sample width and to 
follow closely the electron density profile. On the other 
hand, there is a simple classical argument supporting 
Chang's conjecture. If we describe the magnetotransport 
in the inhomogeneous Hall bar by a local Ohm's law with 
a position-dependent resistivity tensor p, we find that the 
current density is largest where the longitudinal resistiv- 
ity is smallest, i.e., along the incompressible strips, if we 
assume that p(y) depends on the local filling factor in 
the same manner as for a homogeneous sample. Such 
a classical local magnetotransport model has previously 
been used to explain, for inhomogeneous samples, a finite 
width of quantum Hall plateaus without the assumption 
of localized states, and to simulate the magnetic-field de- 
pendence of Hall-type voltages in Hall bars with internal 
contacts. 19 ' 20 More recently such a local model has also 
proven useful for the understanding of the current and 
electric field distribution in antidot systems close to the 
breakdown of the QHE. 21 ' 8 

Since the dissipative non-equilibrium current will 
lead to a position-dependent electrochemical poten- 
tial, which changes the equilibrium electron-density and 
electrostatic-potential distribution, we have to generalize 
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the Thomas- Fermi-Poisson scheme for the self-consistent 
calculation of the latter to include the current-induced 
changes. This will be done in sect. II. Typical results 
obtained from this approach are presented in sect. III. 
In sect. IV we will discuss the most relevant results in 
the light of the motivating experiments 10 , and we will 
comment on apparent limitations of our model. Some 
preliminary results of this work have been presented 
previously. 22 



II. MODEL 

Following previous work 14 ' 18 we model the Hall bar 
as a 2DEG in the plane z = 0, being restricted to the 
strip \x\ < d and translation invariant in the y direction. 
Physically, the confinement potential of the 2DEG is pro- 
duced by a homogeneous background charge in the strip, 
so that the density of free charges has the form p(x)8(z) 
with 



p(x) — e[n — n e i(x)] 0(d 2 — x 2 



(1) 



with n e \(x) and no the surface densities of 2DEG and 
background, respectively. Electrostatic boundary condi- 
tions are fixed by the assumption of metallic halfplancs 
of constant potentials Vl and Vr in z = 0, x < —d, and 
in z = 0, x > d, respectively, and by dielectric constants 
k> and n < in the half spaces z > and z < 0. With 
these assumptions, for z ^ the electrostatic potential 
V(x, z) = Im_F(£) satisfies the Laplace equation and can 
be written as the imaginary part of an analytic function 
F(() of the complex variable ( — x + iz, which is deter- 
mined by its boundary conditions on the real axis. 23 ' 13 ' 14 
We shall recall briefly the use of this electrostatic 
model for the description of the thermodynamic equilib- 
rium state 15 ' 18 and then propose an extension to current- 
carrying stationary non-equilibrium states. 



A. Thermal equilibrium 
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(4) 



where w(Q = ^/ d 2 — £ 2 is analytic except on cuts at 
z = 0, \x\ > d, and R = (k> + k<)/2. By integration one 
obtains the (bare) confinement potential for n e \(x) = 
as 



V a (x) = -E Q ^jl-{^y , Eo = ^-n d. (5) 



The Hartree potential due to the 2DEG follows as 

r-l 
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with the kernel 



In 



(6) 



(7) 



The electron density is, in turn, determined by the po- 
tential V(x). We assume that V(x) varies slowly on the 
scale of typical quantum mechanical lengths, notably the 
magnetic length l m — yjfi/mu> c defined by the cyclotron 
frequency ui c = eB/mc, and calculate the electron den- 
sity in the Thomas-Fermi approximation 



n cl (x)= J dED(E)f(E + V(x)-p,*), 



(8) 



where f(E) = 1/[1 + cxp(^/fcsT)] is the Fermi distri- 
bution and D{E) the (single particle) density of states 
(DOS) of the 2DEG. For the Hall bar in the absence of 
a magnetic field, B — 0, we take D(E) = D 0(E), with 
Do = m/tirh 2 ), while for large B we will use a (suitably 
broadened) Landau DOS. 

Solving Eqs. (2) and (8) self-consistently for constant 
electrochemical potential p* , we obtain the electron den- 
sity and the electrostatic potential in the thermal equi- 
librium state. The value of p* determines the average 
electron density and vice versa (for fixed T, B, etc.). 



To determine n e i (x) we need the electrostatic potential 
V(x) = V(x, z = 0) in \x\ < d. It can be written as the 
sum of three terms, 



V(x)=V (x) + V e (x) + V H (x), 



where 



vl , Vl + Vr , Vr-Vl 
VAx) = 1 arcsin , 

SW 2 7T \d 



(2) 



(3) 



is determined by the potential values on the in-plane 
gates and accompanied by compensating induced charges 
on the gates. The free charges determine V(x, z) = 
ImF(C) with 18 [we write the potential as potential en- 
ergy of an electron, i.e., including a factor — e] 



B. Local equilibrium with imposed current 

If a stationary net current is imposed on the Hall bar, 
this leads to position-dependent current densities and 
electric fields which, in the linear response regime, are 
interrelated by Ohm's law. The relevant field driving 
the net current is the gradient of the electrochemical 
potential, Vp*, which vanishes in thermal equilibrium. 
In view of Eq. (8), we should expect that a position- 
dependent p* will lead to a modified electron density 
n e i(x) and this, according to Poisson's equation, to a 
modified electrostatic potential. Our aim is to extend 
the self-consistent Thomas-Fermi-Poisson approximation 
for the thermal equilibrium state to the stationary non- 
equilibrium situation defined by the imposed current. To 
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do so, we adopt the widely used assumption of local equi- 
librium: thermodynamic variables are assumed to vary 
only slowly in space and to satisfy locally the same rela- 
tions as they would do in a homogeneous thermodynamic 
equilibrium state. Then, for a given position-dependent 
p*(x), we can use again Eqs. (2) and (8) to calculate 
V(x) and n e i(x). 

We will now formulate a model that allows to calcu- 
late V/U*, and thus the position-dependent electrochemi- 
cal potential up to a constant, provided the electron den- 
sity n e i(x) and the total current 



dxj y (x,y) 



(9) 



are given. In the spirit of the Thomas- Fermi and the local 
equilibrium approximation, which assume that electro- 
static potentials and thermodynamic variables vary on a 
scale much larger than quantum lengths, we also assume 
that the current density j(r) — (jx(x,y),j y (x,y)) and 
the electric field satisfy the local version of Ohm's law, 



p(r)j(r) =E(r) = Vp*(r)/e, 



(10) 



where the resistance tensor p(r) = [<7(n c i(r))] 1 is as- 
sumed to depend on position only via the electron density 
riei(r). 

Assuming a stationary situation with translation in- 
variance in y direction (i.e. p, j, and E are independent 
of y), we obtain from the Maxwell equations V • j(r) =0 
and VxE(r) = that the current density j x across and 
the field E y along the bar are independent of x, 



With Pyy = p x 

further get 

3y{x) 



j x = 0, E y (x)=E° y . (11) 
pi{x), p xy = -p yx = ph(x) we 



Pl{x 

From these results we obtain 



El Ex (x) = ^E°. (12) 



V*{x,y) = Po+eEl 



dx' 



Pl{x) 
, Ph(x') 



Pi{x') 



(13) 



where p^ occurs as an undefined constant, and from 
Eq. (9) we get the normalization 



J_ d dx[l/pi(x)} 



(14) 



The new self-consistency problem for the stationary 
state is completely defined, if we choose a model for the 
dependence of the conductivity tensor on the electron 
density. To solve it iteratively, we start with zero current, 
1 = 0, and solve the old equilibrium problem. Then we 
take a fixed value / ^ and proceed as follows. 

In the next step of the new iteration, we use n e \{x) 
from the previous step and calculate for the given / the 



electrochemical potential (13). Then we put this into 
Eq. (8). To compensate the y- linear term, we add an 
identical term to the electrostatic potential. This guar- 
antees that the electron density remains independent of 
y. Then we solve the "old" problem, Eqs. (2) and (8) 
with the modified x-dependent electrochemical potential, 
self-consistently to determine V(x) and n c \{x), choosing 
Pq so that the average electron density remains the same 
as without current. Convergence of the "old" problem 
completes this step of the new iteration. The steps of 
the new iteration are repeated until there is practically 
no further change of n e \(x). 

We have performed these self-consistent calculations 
for two types of electrostatic boundary conditions. First 
we started with a symmetric electron profile and Vq = 
Vr — Vl = 0, and assumed that an applied current does 
not change the potentials of the in-plane gates. Then, in 
a strong magnetic field and under an imposed current, 
we obtain an asymmetric electron profile that is shifted 
towards one of the sample edges. This is an obvious 
consequence of the Lorentz force on the drifting electrons. 
In this situation, the applied current leads to a change of 
the induced charges in the in-plane gates. 

Since dV(x, z)/dz = Re{dF/d(], we can use Eq. (4) 
to calculate the induced charge density in the in-plane 
gates, pind(^) = [K/(27re)]ReF'(x + «0 + ). For the total 
induced charge Qr = J d dxp- m d{x) in the right gate this 
yields 



Qr 
en d 



d£ 



n e \(^d) 



n 



arctan J j — - . (15) 



The corresponding result Q l for the left gate is obtained 
from Eq. (15) by replacing £ under the square root by — £. 
The sum of these induced charges compensates the free 
charges, Qr + Ql = — dxp(x), and their difference 
vanishes if the electron density is symmetric, n e i(— x) = 
n e \(x). 

The asymmetric density profile, resulting for the 
current-carrying stationary state from the requirement 
Vq = 0, leads to Qr ^ Ql- Since Qr + Ql is kept con- 
stant, this means that charge must flow from one gate to 
the other as the stationary state is established. As an al- 
ternative electrostatic boundary condition, which may be 
more realistic in certain situations, we investigated the 
"floating gate" condition, requiring that Qr (and thus 
Ql) is kept constant and a finite voltage between the 
in-plane gates builds up. 



III. RESULTS 

A. Classical regime 

In the classical regime of low magnetic field and high 
temperature, the magnetic field should not affect the 
thermodynamic equilibrium state 24 (Bohr-van-Leeuwen 
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theorem). Therefore, we use the simple DOS D(E) = 
Do 0(E), which for zero temperature renders Eq. (6) into 



Vff(£d) 



d('K(t,0[E F -V(C'd)], (16) 



where the dimensionless parameters (3l < Pr define the 
edges of the density profile, V(@Ld) = V(0 R d) = Ep = 
pf(T = 0), and where 



a = Trn / (E D ) = irao/d , 



(17) 



with ao = Rfi 2 / (2me 2 ) the screening length. [Due to a 
misprint, in the denominator in Eq. (15) of Ref. 18 the 
factor 7r is missing.] Together with Eq. (2), Eq. (16) rep- 
resents a linear integral equation for V(t;d) in the interval 
/3l < £ < (3r, which can easily be solved numerically. So- 
lutions of this linear equation are used as starting points 
for all numerical calculations at finite temperature and 
magnetic field, which lead to non-linear integral equa- 
tions that must be solved iteratively. 

As has been shown in Ref. 18, with decreasing values of 
a screening becomes more effective, and a voltage Vq ap- 
plied across the bar leads to a shift and deformation of the 
equilibrium electron density profile n c \(x). At finite tem- 
perature, the sharp edges of the zero-temperature profiles 
are smeared out. 

In the non-equilibrium calculations, we use the Drude 
model for the resistivity tensor, 



pi(x) = l/a (x), ph{x) = lo c t pi(x), 



(18) 



with <jq(x) — (e 2 r /m)n B \(x). Then, according to 
Eq. (12), the current density j y (x) = <r (x)Ey is pro- 
portional to the electron density, and E x (x) = lu c t E® 
is constant. Strictly speaking the last result, which 
follows from the fundamental linear response equation 
jx(x) — a xx (x) [E x — uj c t E y ] = 0, can be justified from 
our local approach only at positions x where n c \(x) ^ 0. 
On the other hand, within our local equilibrium model 
we expect the electrochemical potential to be constant in 
regions where no (dissipative) current flows and where no 
electrons are. Therefore we put E x (x) = if n e \(x) = 0. 
In practical calculations we use 



E x (x) = lo c t E^0(n c i(x)/n Q - e) 



(19) 



with e = 10~ 4 , which defines effective edges < bR of 
the density profile by n i(&i) = n c \(bn) = en . 

Typical results are shown in Fig. 1, where for the 
current-carrying states the variation of the electrochemi- 
cal potential, Ap* = p*(b R ) - p*{b L ) = (b R - b L )uj c TE°, 
is fixed as uh — A/i* /Eq = 0.2. This yields the current 
I = uh [2d/(&a — &l)] (e 2 fi e i/mu> c ) Eq / e, where 2dn e \ = 
f_ d dxn c i(x) is fixed, n c i/n = 0.6278. We observed 
that, for fixed vanishing gate voltage and < uh 0.2 
the shift (and deformation) of the density profile in- 
creases roughly linearly with uh- Under floating gate 



conditions, for uh = 0.2 a voltage Vr — Vl — O.SlOE'o 
builds up between the in-plane gates. This is larger than 
the linear extrapolation uh [2d/(bR — i>z,)]l?o due to the 
singular slopes of the self-consistent potential at the edges 
of the gates (see Fig. 1). Figure 1 shows that, over the 
whole range of finite density, the self-consistently calcu- 
lated total potential V(x) follows very closely the lin- 
ear position-dependence of the electrochemical potential. 
This observation holds for both the fixed-gate-voltage 
and the floating-gate boundary conditions. In the follow- 
ing we will consider only the more realistic floating-gate 
boundary conditions. 
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FIG. 1. a) Density profile n B \(x) and b) total potential 
V(x) in equilibrium (gray lines) and with current I oc uh 
(black lines) for fixed gate voltage Vr = —Vl — (solid) 
and for floating gates (dash-dotted lines), and c) induced 
Hartree potentials V(x; I) — V(x, 0) for both cases. The thin 
lines indicate the electrochemical potentials p*(x). a = 0.01, 
#r = -Pl = 0.848. 



B. Quantum regime 

For strong, quantizing magnetic fields we should use a 
suitable form of a broadened Landau DOS 



D{E) 



2nll 



J2 A ^(E), 



(20) 



n=0 



where A n (E) is the spectral function of the n-th Lan- 
dau level with energy eigenvalue e n — Tiuj c (n + 1/2), and 
g s = 2 accounts for spin degeneracy. The model to be 
used for the resistivity tensor should show the charac- 
teristic behavior known from the quantum Hall regime, 
notably (nearly) vanishing pi(x) at (even) integer local 
filling factor v{x) = 2 r Kl m n e \(x). Moreover, the approxi- 
mations for the conductivity tensor and the DOS should 
satisfy certain consistency relations (consequences of the 
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equation of continuity). 27 However, before we address 
such sophisticated questions we want to present a very 
simple model. 



1. Simplified model 

First we follow previous work 15,18 and consider the 
bare Landau DOS, taking in Eq. (20) A n (E) = 5{e n -E). 
This leads, in thermal equilibrium, to the appearance of 
incompressible strips of finite width at integer values of 
the local filling factor. The temperature dependence of 
these strips has been discussed in Ref. 15 and the depen- 
dence of position and width of the strips on magnetic field 
and applied gate voltage Vr — Vl has been investigated 
in Refs. 18 and 25 for the present Hall bar geometry. To 
investigate the effect of an imposed current, we first use a 
simplistic model for the resistance tensor, that has been 
used successfully for the calculation of the current den- 
sity in an antidot system in a strong magnetic field. 26 ' 21 
For <jy X (x) = —o xy (x) = (th{x) we take 



a H (x) = (e 2 /h) v{x) , 



(21) 



which yields the correct values at integer filling factors, 
but no quantum Hall plateaus. To simulate the behav- 



oi near integer filling v 



ior of a xx 

approximate 21 the longitudinal conductivity as 

o-i(x)=o- H (x)[e+(2~v(x)) 2 /4] : (22) 



-31 

screened potential 
_32 L [ electro- 
Landau I 




FIG. 2. a) Self-consistent potential V(x) and b) normalized 
density v(x) for zero (thick solid gray) and finite (thick solid 
black lines) current, calculated for model (22) with e = 0.002. 
a) also shows the corresponding n*{x) (dash-dotted) and the 
Landau levels V(x) + Q(n + 1/2) (dashed lines), and b) 
the current density (long-dashed line, arbitrary units), c) 
shows the current-induced change of the self-consistent po- 
tential (black). Model parameters (see text) are a = 0.02, 
Q = hui c = So/200, k B T/Q, = 0.04, U H = A^* = 3«. 



with a small but finite positive value e (~ 10~ 3 ) to avoid 
divergencies. This describes correctly that pi(x) oc <?i(x) 
becomes very small at the incompressible strips with local 
filling factor v(x) — 2, although the analytical dependen- 
cies for v ^ 2 are not correct (see Fig. 3b below). 

Nevertheless, this simple model is able to reproduce 
characteristic features observed in the experiment 10 as 
is shown in Fig. 2. The current flows preferably along 
the incompressible strips, where v(x) = 2 and the lon- 
gitudinal resistivity is smallest, and there the gradient 
of fj,*(x) is largest. Due to the Thomas-Fermi- Poisson 
self-consistency requirement, the total potential V{x) is 
forced to follow /it* (x) closely, so that the current-induced 
change of the electron density profile is small (which 
keeps the change of electrostatic energy small). The in- 
duced AV(x) follows closely n*{x) and varies mainly in 
the region of incompressible strips. 



2. Gaussian level broadening 

As a more realistic model for the longitudinal conduc- 
tivity we use the Gaussian model 27 



e 2 g s 



dE 



dE 



n=0 

with the spectral function 



OO -. 

J2( n +2)[^ T ME)] 2 (23) 



A n {E) 



exp( 



£;] 2 /r 2 



(24) 




FIG. 3. Longitudinal conductivity 07 versus filling factor v 
for the Gaussian model (23) with F/huc = 0.035 for (a) high 
and (b) low values of the reduced temperature t — ksT '/hui c . 
The thick solid curves coincide with the limits for high tem- 
perature [fcflT > Huj c , (a)] and zero temperature [fcsT <C T, 
(b)]. The long-dashed curve in (b) indicates the model (22). 
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which then, for consistency reasons, should also be used 
in the DOS, Eq. (20). An alternative, which leads to 
qualitatively the same results, would be to use the self- 
consistent Born approximation 28,20 which would replace 
the normalized Gaussians by normalized half-ellipses. To 
avoid divergencies, we replace ai of Eq. (23) for v > 1 by 
max((T(, C///10 4 ). 

Together with v = 2^ J dE f(E- fj,)D(E), Eqs. (20), 
(23), and (24) can be used to calculate ai as function of v. 
Results are plotted in Fig. 3 for several temperatures. For 
high temperatures, fcgT > 0.3 fko c , one gets the modified 
Drude result, 07 = <jh/ (w c T gauss ) with an = (e 2 /h) v and 

ft/ r gauss = \J 7r/2r. 

To proceed, we first investigate the effect of the 
Landau-level broadening T on the existence and width 
of the incompressible strips. Results of self-consistent 
calculations are shown in Fig. 4. The width of the incom- 
pressible strips shrinks with increasing temperature and 
with increasing level broadening. For ksT '/hco c < 0.04 
and T/fujj c < 0.1 clearly visible incompressible strips ex- 
ist. Thus, collision broadening of the Landau DOS does 
not change the screening properties of the 2DEG quali- 
tatively, provided the width of the Landau levels remains 
small enough as compared with the cyclotron energy. 
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FIG. 4. Density profile [a), b)] and potential [c), d)] cal- 
culated with a Gaussian DOS; b) and d) show results in the 
region of the incompressible strip for several values of tem- 
perature, ksT /hid c = t/50, and level broadening, 7 = Y /Tiuj c . 
The thick curves are plotted in a) and c) for one half of the 
symmetric sample (a = 0.02, hu c /Eo = 0.005). 

Next we perform the self-consistent calculation of the 
charge and current densities and of the electrostatic and 
electrochemical potentials for the Gaussian model. To 
achieve convergence of the nested self-consistency loops 
for given values of temperature T, cyclotron energy f2 = 
huj c and total current / cx Ur = Afi* , we proceed as fol- 
lows. First we define the density profile by solving the lin- 
ear integral equation [Eqs. (2) and (8)] for T = 0, B = 



and I — 0. Then we raise, still for B = and 1 = 0, the 
temperature stepwise up to the value fc^T = 0.3f2 and 
solve at each step the non-linear problem iteratively using 
a Newton-Raphson procedure. At this high temperature 
all quantum effects are smeared out, and we can replace 
the B = DOS by the Gaussian Landau DOS corre- 
sponding to the required Q value without convergence 
problems. Now we raise stepwise the current until the 
required value is reached. This calculation is equivalent 
to the solution of the Drude problem discussed above. 
When self-consistency is achieved, we lower the tempera- 
ture stepwise until the required (low) value is reached. In 
each step we iterate until full self-consistency is achieved, 
using the previous potential profile and the conductivity 
tensor with the density profile of the previous step as 
starting conditions. 
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FIG. 5. Self-consistent results for a) filling factor v(x), 
b) current density j y (x), and c) electrostatic and elec- 
trochemical potentials, V(x) and n*(x), at five tempera- 
tures t = ksT/liuic, calculated from the Gaussian model, 
(a = 0.02, huj c /E = 0.0053, T/hoj c = 0.03, U H = fiw c ). 

Figure 5 shows the self-consistent results for several 
intermediate temperatures. At the highest temperature 
(dashed lines) one observes Drudc-like behavior: the cur- 
rent density [Fig. 5(b)] is proportional to the electron 
density [note that Fig. 5(a) shows the latter only near lo- 
cal filling factor v(x) = 2, while i/(0) = 2.25] and the clec- 
trostatical and electrochemical potential increase nearly 
linearly across the 2DEG. With decreasing temperature 
the 2DEG develops incompressible strips with low longi- 
tudinal resistivity and the current density is increasingly 
confined to the incompressible regions. Simultaneously 
the potentials develop a steplike behavior with variation 
across the incompressible strips and plateaus in the com- 
pressible regions. 

To evaluate the current-induced electrostatic poten- 
tial AV, we perform the self-consistent calculation 



G 



with and without applied current and define AV(x) = 
V(x; I, B, T) - V(x; 0, B, T). A typical result is shown 
in Fig. 6. The main difference between this result and 
Fig. 2 is that now the current density is confined more 
strictly to a narrow region along the incompressible strips 
(see dashed lines in the middle panels of the figures). 
The more rapid decrease of the current density from the 
large values in the incompressible strips to the small val- 
ues in the compressible regions is caused mainly by the 
much steeper increase of <Ji(v) with increasing \ v — 2| [see 
Fig. 3(b)]. The ratio between the values of the current 
density in the compressible region and those in the in- 
compressible strips is also smaller, since we used near 
i/ = 2a smaller cutoff e = min[<7;(i/)/(7fl-(z/)] in the Gaus- 
sian model (e = 1(T 4 ) than in model (22) (e = 2x 1(T 3 ). 
As a consequence, the variation of the electrochemical 
potential (dash-dotted line of upper panel) and of the 
current-induced electrostatic potential (lower panel) is 
practically confined to the region of the incompressible 
strips. 
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FIG. 6. As Fig. 2, but with oi calculated from the Gaus- 
sian model (23) instead of Eq. (22), for U H = Ap* = fl 
(r/huj c — 0.03, all other parameters as in Fig. 2). 

The width and position of the incompressible strips 
and thus the locations of strong variation of the current- 
induced potential change strongly with varying magnetic 
field, i.e., with varying filling factors of the Landau levels. 
In Fig. 7 we show results for selected values of the mag- 
netic field, leading to filling factors in the center of the 
Hall bar that vary between z/(0) = 1.62 and i/(0) = 4.52. 
The temperature is always chosen so low that the incom- 
pressible strips are well developed (fcBT/?iw c < 0.04). 

These results are easily understood. For Huj c > 5.7 x 
10~ 3 Eq no incompressible strips exist, v(x) < 2 for all 
|x| < d, and the current density is largest near the center 
of the sample, where the filling factor is largest and the 
longitudinal resistivity (pi oc <ji) is smallest. Therefore 



the gradient of AV(x) is largest in the center of the Hall 
bar. If i/(0) < 2 is very close to 2, pi(x = 0) is very 
small, the current density has a sharp maximum in the 
center, and the potential profile has a strongly non-linear 
appearance ("type II" behavior). If i/(0) becomes consid- 
erably smaller than 2, pi(x) has a broad maximum near 
x = and the current density profile follows essentially 
the density profile, similar to the Drude case. This leads 
to an essentially linear potential profile ("type I"), as is 
seen in the top curve of Fig. 7 for Uui c = 7.0 x 10 _3 -E , 
with i/(0) = 1.62. 
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FIG. 7. Current-induced part AV(x) of the self-consistent- 
ly calculated electrostatic potential in units of the cyclotron 
energy Q — hui c , for several values of f2. The numbers in 
the figure indicate (Q/Eo) x 10 3 , those on the right hand side 
the corresponding values of f(0). For clarity, the traces are 
shifted vertically by an arbitrary amount. 

For slightly lower magnetic field, an incompressible 
strip with z/(0) = 2 and (nearly) vanishing pi occurs in 
the center. Then the variation of AV[x) is confined es- 
sentially to this strip (type II) . With still decreasing Uuj c 
this strip splits into two, which move with decreasing 
magnetic field towards to edges of the Hall bar. The elec- 
trochemical potential and (apart from some minor edge 
effects) the current-induced electrostatic potential then 
drop only across these incompressible strips ( "type III" ) , 
as is seen in Fig. 7 for the curves with 2.04 < i/(0) < 3.77. 
As with further decreasing magnetic field the filling factor 
in the center region comes close to 4 and thus pi(x = 0) 
becomes small, a considerable part of the current flows 
through this center region. Since at the same time the 
strips with local filling factor 2 become very narrow, for 
i/(0) ^ 4 a considerable part of the induced potential 
drops in a broad center region (type II). For i/(0) > 4 the 
center region becomes again compressible, with constant 
AV(x), and the incompressible strips with v(x) = 4, 
across which now most of the Hall voltage drops, move 
away from the center (type III). The lowest trace in Fig. 7 
shows such a situation with tiny structures at the edges 
of the electron density profile which are due to the outer 
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incompressible strips with local filling factor 2. 

IV. CONCLUSIONS 

Our results for the current-induced Hall-potential pro- 
file (Fig. 7) reproduce characteristic features of the ex- 
periment of Ahlswede et al. 10 . If the filling factors j^(O) in 
the center of the sample are close to, but below integer 
values, the potential drops in a non-linear fashion in a 
broad center region. For v(0) values slightly larger than 
the integer values, the potential is constant in the center 
region and drops exclusively across the incompressible 
strips. Of course, our results show this characteristic be- 
havior only near even integer values of ^(0), since we have 
neglected spin-splitting, whereas in the experiment spin- 
splitting is resolved and this behavior occurs also near 
small odd- integer values of ^(0). 

This characteristic dependence of the Hall-potential 
profile on the magnetic field cannot be explained by 
the previous calculations assuming dissipationless Hall 
currents 17,18 ' 25 and emphasizes the importance of dissi- 
pation. From the nice qualitative agreement of our re- 
sults with the experimental data we conclude that our 
local equilibrium approach, which combines dissipative 
transport with screening effects and allows to calculate 
electron and current density as well as electrostatic and 
electrochemical potential self-consistently, contains most 
of the relevant physics. There is, however, room and need 
for improvements. 

One desirable improvement concerns the effectiveness 
of narrow incompressible strips. To avoid numerical di- 
vergencies, we used a cut-off e • <jh for er; at even integer 
values of v, with e = 10~ 4 . If we take the limit e — > 
and sufficiently low temperatures, pi{x) becomes expo- 
nentially small in the incompressible strips. Then, ac- 
cording to Eq. (14), the electric field along the Hall bar 
and, therefore, the longitudinal resistance become expo- 
nentially small whenever an incompressible strip exists, 
and not only for a limited interval of magnetic field val- 
ues in a plateau region of the QHE. To eliminate this 
unreasonable behavior, we should include a mechanism 
(other than the simple cutoff) that limits the current den- 
sity in, and thereby the voltage drop across, narrow in- 
compressible strips. Such a mechanism could make the 
incompressible strips with local filling factor v{x) = 2 
ineffective for magnetic fields with i/(0) > 3. This would 
turn the "type III" curves with v(Q) = 3.23 and 3.77 in 
Fig. 7 into quasi-linear "type I" curves and would elim- 
inate the tiny edge-near structures in the three lowest 
curves. Both changes would improve the agreement with 
the experiment. 

Several physical effects may lead to such a mechanism. 
One is the non-local relation between the current density 
and the driving electric field, which we have approxi- 
mated by a strictly local one. Another one is Joule heat- 
ing, which is most effective where the current density 



is high and may destroy narrow incompressible strips, 
i.e., lead to a local breakdown of the QHE. A systematic 
treatment of heating effects will require the considera- 
tion of energy balance and heat conduction, as has re- 
cently been pointed out by Akera 29 in his hydrodynamic 
approach to quantum Hall systems in the breakdown 
regime. A consideration of the heating processes relevant 
under QHE conditions (e.g. "quasi-elastic inter-Landau- 
level scattering") ' may even demand a treatment be- 
yond a local hydrodynamic approximation and require 
a more microscopic non-local description of stationary 
current-carrying non-equilibrium states intermediate be- 
tween the zero-resistance quantum Hall state and the 
finite-resistance breakdown state. 32 Such a microscopic 
approach to heating and resistive processes may also open 
the possibility of a unified description of dissipative cur- 
rents, which we have considered phenomenologically in 
the present paper, and non-dissipative equilibrium cur- 
rents, which we have mentioned in the introduction but 
completely neglected in the calculations. 

Finally we want to mention that in our model calcu- 
lations the imposed current leads to a broadening of the 
incompressible strips on one side of the sample and to 
a narrowing of the corresponding strips on the opposite 
side (of course the strips exchange their role if we invert 
the direction of the current). This asymmetry is clearly 
seen in Fig. 7, and can also be observed in Fig. 1 of 
Rcf. 10. A systematic investigation of this effect may be 
of interest. 
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